Zitterbewegung-like effect near the Dirac point in metamaterials and photonic crystals 
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We present a physical explanation of Zitterbewegung-like effect near the zero-refractive-index point in a 
metamaterial slab in this paper. Between the negative and positive refractive index regions centered at the 
zero-refractive-index point, the transmittance spectrum distribution of the metamaterial slab is asymmetrical. 
When a symmetrical pulse propagates through the metamaterial slab, its transmitted spectrum becomes asym- 
metrical due to the asymmetry of the transmittance spectrum of the slab, leading to a transmitted pulse with 
an asymmetrical temporal shape. The asymmetry manifests a kind of temporally tailed oscillations, i.e., the 
Zitterbewegung-like effect. Further, the effect of the temporal and spatial widths of pulse, and the thickness of 
metamaterial slab on the tailed oscillations of the transmitted pulse has also been discussed. Our results agree 
well with what the other researchers obtained on the strength of relativistic quantum concepts; however, the 
viewpoint of our analysis is classical and irrelevant to relativistic quantum mechanics. 

PACS numbers: 78.20.Ci, 41.20Jb 



I. INTRODUCTION 

The propagation of electromagnetic waves near the Dirac 
point (DP) in photonic crystals and metamaterials has been 
investigated intensively in recent years llll- Ulll . With some 
special structure designs of both kinds of artificial materials, 
two energy bands touch each other, forming a pair of cones, 
namely, so-called "Dirac point" in optics. For a monochro- 
matic electromagnetic wave propagating near the DPs, some 
abnormal transmission properties have been discovered, such 
as pseudo-diffusive scaling and conical diffraction 10]. 

The propagation of optical pulses with the central frequency at 
the DP in a photonic crystal or metamaterial slab has also been 
investigated by some research groups who show that tempo- 
rally tailed oscillations of the transmitted pulses should appear 
llHE]]. 1° meir papers, this kind of oscillation is called Zitter- 
bewegung (ZB) effect for photon. 

The existence of ZB effect was first proposed by Erwin 
Schrodinger for relativistic electrons in free space, and the in- 
terference between the positive and negative energy states was 
considered as the origin lfl2ll . ZB effects in crystals, supercon- 
ductors, semiconductor nanostructures, graphene, and other 
physical systems have also been demonstrated II 13142711 . And 
ZB effect for photon has been proposed in two-dimensional 
photonic crystal for the first time by simulation analysis by 
Zhang 1 8], in which he used a plane wave with Guassian 
shape in time domain, passing through an open slit, to ex- 
cite the photonic crystal slab. In another paper for the case 
of sonic crystals, Zhang and his co-worker demonstrated the 
ZB effect in experiment and mentioned that the origin of 
this effect was the interference between two linear modes 
around the DP [9]. Shortly thereafter, DP in the negative- 
zero-positive index metamaterial (NZPIM) was first proposed 
by Wang et al. II 1 Oil . Thereupon the ZB effect of optical pulses 
in NZPIM was characterized [11]. Both the research groups 
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named the temporally tailed oscillations of the transmitted 
pulses as ZB effect and thus explained it on the strength of 
some concepts of the relativistic quantum mechanics, such as 
DP and Dirac equation. This is undoubtedly an interesting and 
helpful understanding of the underlying mechanism of the ef- 
fects in photonic crystal and metamaterial. However, we think 
that it may be not suitable to use ZB directly as the name of 
the tailed oscillations of optical pulses. So, here, we introduce 
Zitterbewegung-like (ZB-like) instead. 

In this work, based on another point of view which is ir- 
relevant to the relativistic quantum mechanics, we present a 
physical explanation of temporally tailed oscillations (or ZB- 
like effect) of the pulses propagating near the zero-refractive- 
index point (called DP in Refs. O [O]) in NZPIM using 
the spatio-temporal filtering analysis. Due to the asymmet- 
rical transmittance spectrum distribution between both sides 
of the zero-refractive-index point, a symmetrical pulse passed 
through a NZPIM slab undergoes an asymmetrical filtering 
effect, which leads to the asymmetrical temporal shape of the 
transmitted pulse, in other words, the tailed oscillations. Fur- 
ther, the effect of the temporal and spatial widths of pulse, and 
the thickness of the NZPIM slab on the tailed oscillations of 
the transmitted pulse can be characterized well. The analysis 
method we use is similar to that of previous works on 
dealing with some detail problems, however, our viewpoint, 
based on classical spectral filtering analysis that is not related 
to the relativistic quantum mechanics, is totally different from 
theirs. Even more interesting is the fact that the results we ob- 
tained are in good agreement with those in Refs. |8J] and 111 111 . 
Using this analysis method, ZB-like effect in two-dimensional 
photonic crystals could also be explained. 



II. MODEL AND THEORY 

We take the homogenous NZPIM slab as the analysis sam- 
ple. This kind of metamaterial has already been demonstrated 
by theoretical and experimental researches based on liquid 
crystals and fi-shaped microstructures from GHz to visible 
regions j28T[32ll . For simplicity, Drude model is chosen as 
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FIG. 1: Schematic of the NZPIM slab in vacuum with thickness h 
along z direction. It is assumed that the slab is infinite along x and y 
direction. 



the parameter of the permittivity and the permeability of the 
NZPIM slab OHH: 
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where a> pe and CO p m 3X6 controllable electronic and magnetic 
plasma frequencies, and y e and y,„ are the damping rates. 



As shown in 13011 . with suitable structure parameters, a> p = 
u) pe = u> pm can be achieved in a back-to-back £2-shaped meta- 
material in GHz region. When a> — a> p = u> pe = cj pm and 
7 - le — Jm ^ Mp, e(a>) and fi(u>) are nearly equal to zero, si- 
multaneously. Thus, the effective refractive index of the sam- 
ple varies from negative, through zero, then to positive with 
the increase of frequency. In such case, the so-called DP is 
the zero point of effective refractive index. 

We consider a paraxial, coherent, TE-polarized, spatially 
and temporally Guassian pulse injecting into the sample along 
z direction in vacuum as shown in the Fig. [T] It is assumed 
that the slab sample is infinite along x and y direction. The 
refractive index is n\(a>) = yj ' e(<x>)ii{<jj) — 1 - lo p /(lo 2 + iya>) 
and the thickness of the sample is h. The function of the pulse 
at the initial plane z — Zo < [while plane z — is the interface 
between vacuum («o — 1) and the NZPIM slab sample] is 
given as 



Ei (x, t, zo) = exp 
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This pulse can be decomposed into its monochromatic Fourier 
components and further be expanded in angular spectra of 
monochromatic plane-wave components 
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where V is the temporal half-width, Wo is the spatial half- 
width, k x = («ow/c)sin#, c is the light velocity in vacuum, 
and (l>o is the central angular frequency of the pulse. 

The transmission coefficient t(8, a>) = \t(9, u>)\ exp[iip(8, ai)] 
of the sample can be deduced from the electromagnetic conti- 
nuity conditions of the two interfaces (plane z — and z = h in 
Fig.[TJ 0411 . In this paper, for the TE-polarized incident wave, 
t(8, (d) can be written as 
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FIG. 2: (Color online) (a) The spatio-temporal spectrum of incident 
Guassian pulse with T = 1.0ns, Wo = IOAq (Aq = 0.03m), and 
f = 10GHz. (b) Distribution of the spectral transmittance (\t(0,a>)\) 
for a sample with h = 40A and y = y e = y m = 10 4 Hz, the blue 
part denotes the stop-band and the rest denotes the pass-band. When 
a> = ai,, = 2/r x 10GHz, the refractive index is zero, (c) The spatio- 
temporal spectrum of the transmitted pulse. 



Here, 0\ = arcsin[«o sin 0/ni(a>)] is the refractive angle of the 
angular spectrum component in the sample which can be ob- 
tained from Snell's law and 6 = w«i(w)/!cos0i/c stands for 
the pha se change when electromagnetic wave cross the sam- 
ple in. 

In this case, the transmission coefficient is equivalent to the 
frequency response of a linear filter 135143911 . Therefore the 
spectral function at the exit end can be written as E,(9, to, z) = 
Ej(8, oj, 0)t(8, to). Ei(6, to, 0) = E0, to, Zo) exp(-//bz ) is the 
function of spatio-temporal spectrum of the pulse arriving at 
the plane z — 0, where k z = ko cos 6 (ko = noa>/c is the wave 
vector in vacuum). We take the spatial and temporal half- 
width of the incident pulse as 10/lo (Aq = c/f is central wave- 
length of the pulse in vacuum) and Ins, respectively, with the 
central frequency f (=10GHz) for an example. The spatio- 
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temporal spectrum of the incident pulse and the transmittance 
spectrum distribution \t{6, (o)\ with slab thickness h = AOAq 
are plotted in Figs. [2 a) and 12b), respectively. Unambigu- 
ously, the output spectrum is tailored by the transmittance 
spectrum as shown in Fig. |2jc). Some spatio-temporal spec- 
trum components are filtered by the stop-band of the transmit- 
tance spectrum. In order to obtain the spatio-temporal shape 
E r (x, t, z) of the transmitted pulse, the inverse Fourier trans- 
form of E,(0, 10, z) has been calculated. But it is difficult to 
solve such a complicated expression analytically, so numeri- 
cal integration method is used here. In our analysis, we only 
need to consider the temporal shape of the transmitted pulses 
[we consider E,(x = 0, t, z)] because ZB-like effect is the tem- 
porally tailed oscillations. 



in. RESULTS AND DISCUSSIONS 

As is shown in Fig.|2jb), the pass-band of the transmittance 
spectrum gradually shrinks to the zero-refractive-index point 
from both sides of it, forming a double-cone band. The un- 
derlying mechanism is that toward the zero-refractive-index 
point from both sides, the critical angle of the total reflection 
at the interface z — decreases, thereby only normal incident 
angular spectra of the pulse can pass through the sample at 
the zero-refractive-index point, which can be deduced from 
Snell's law. Even more importantly, the transmittance spec- 
trum is asymmetrical between negative and positive refractive 
index regions in respect that ni(a>) — 1 - upjjiix) 2 + iyai) is 
asymmetrical for the two regions. The asymmetrical trans- 
mittance spectrum gives an asymmetrical filtering effect to the 
incident spatio-temporal spectrum. This asymmetry in fre- 
quency domain is the origin of the tailed oscillations in time 
domain. In a sense, a pulse shape with tailed oscillations could 
be seen as a kind of asymmetry. At the same time, the trans- 
mittance spectrum is symmetrical for negative and positive in- 
cident angle 9, deduced from Eq. (|4]i. Thus, the spatial shape 
of the pulse is symmetrical (seen Fig.|2]in Ref. lUllO . In ad- 
dition, it is worth noting that the temporally tailed oscillations 
will also arise when the central frequency of the pulse diverge 
from the zero-refractive-index point slightly, on condition that 
the asymmetrically filtering effect still holds. 

Equation (|4]i determines the sample's transmittance spec- 
trum that is a function of the sample thickness h. As h in- 
creases, the transition region from the pass-band to stop-band 
becomes steeper and steeper with increasingly severe oscilla- 
tions, which will influence the filtering effect and, then, the 
oscillating properties of the transmitted pulse. In Fig.|3a), we 
plot the relative intensity \(E,(x = Q,t,z)/E,(x = Q,t,z) m ax\ 2 
of the transmitted pulses with the sample thickness h = lO/lo, 
30/lo, 50/lo, 70 Ao, and 90/lo, respectively, while F = 1 .0ns and 
Wo = lO/lo are fixed. The initial center of the pulse is assumed 
at zq - —20Aq. The oscillating strength (the ratios of relative 
intensities between the second peak and the first peak) of the 
tailed oscillations increases with h while the oscillating period 
T (time spacing between the second peak and the first peak) 
remains nearly unchanged, which can be seen from Figs. |3jb) 
and[3]c). When h increases to several hundreds of Aq, the dis- 
persion effect is not negligible and should seriously distort the 
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FIG. 3: (Color online) (a) Relative intensity of the transmitted pulse 
from the calculating data of the inverse Fourier transform with sam- 
ple thickness 10/lo, 30 >Aq, 50 'Ao, 70 Ao, and 90 Ao, respectively, with 
the fixed T = 1.0ns and W = 10A . (b) Oscillating strength versus h. 
(c) Oscillating frequency l/T versus h. 



temporal shape of the pulse. 

Pulse with different temporal half-width F will undergo dif- 
ferent responses when propagating in the NZPIM slab sam- 
ple. The incident spatio-temporal spectrum varies with T, thus 
results in different filtering effects. Near the zero-refractive- 
index point, the dispersion is considered to be approximately 
linear lHoll . however, the linear dispersion approximation is 
not valid when T becomes very small (few optical cycles in 
time domain but extremely large bandwidth in frequency do- 
main). In Fig. Ufa), we plot the relative intensity of the trans- 
mitted pulses with F = 0.5ns, 1.0ns, and 1.5ns, respectively, 
while Wq = 10 Ao and h - AOAq are fixed. Figure |4|b) il- 
lustrates the relationship between the oscillating strength and 
T. As T < 0.8ns (only several optical cycles), the oscillating 
strength varies sharply due to the influence of high-order dis- 
persions. While for F > 0.8ns, the high-order dispersions are 
negligible, so the oscillating strength has a relatively flat de- 
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FIG. 4: (Color online) (a) Relative intensity of the transmitted pulse 
with different temporal half-width 0.5ns, 1.0ns, and 1.5ns, respec- 
tively, with Wo = 10/io and h = 40 Aq. (b) Oscillating strength versus 
T. (c) Oscillating frequency 1/T versus 1/T. 



FIG. 5: (Color online) (a) Relative intensity of the transmitted pulse 
with different spatial half- width 10/lo, 20/to, and 30/lo, respectively, 
with T = Ins and h = 40 Aq. (b) Oscillating strength versus Wo. (c) 
Oscillating frequency 1/T versus 1/Wq. 



pendence on F in the framework of filtering effect alone. The 
oscillating frequency 1/T increases linearly with the 1/T un- 
der the linear dispersion approximation as shown in Fig.[4fc) 
(when 1/T > 1.25GHz, the dispersion is no longer linear). 

Though the ZB-like effect occurs in time domain, the spa- 
tial width of the pulse does, indeed, influence the oscillating 
properties. When Wo increases, the spatial part of the incident 
spatio-temporal spectrum becomes narrower and narrower (6- 
axis direction in Fig. [2 a)). If we consider the limiting situa- 
tion, that is, Wo tends to infinity, the pulse only has the angular 
spectra in normal incident direction, under the circumstances, 
no spectrum component should be filtered by the transmit- 
tance spectrum. So the temporally tailed oscillations would be 
too weak (nearly vanished) to observe as long as Wo becomes 
enough large. In Fig. |3a), we plot the relative intensity of 
the transmitted pulses with Wo = 10 Aq, 30 Aq, and 50 Aq, re- 
spectively, while r = Ins and h = 40/lo are fixed. Seen from 
Fig-Bb), we know the oscillating strength decreases with the 



W). The oscillating frequency 1/T raises sharply with small 
values of 1/Wo and more flatly after coming up to a certain 
value (about 0.06, see Fig.|5jc)). 

It should be pointed out that, although we have only dis- 
cussed the optical ZB-like effect in metamaterial in the above 
analysis, the similar phenomena in two-dimensional photonic 
crystals could be interpreted in the same way. In a recent pub- 
lication 1 4], an extremely similar transmittance spectrum in 
two-dimensional hexagonal photonic crystals near the DP has 
been revealed. 



IV. CONCLUSIONS 

In conclusion, based on the spatio-temporal filtering analy- 
sis, we have presented a physical explanation for the ZB-like 
effect in NZPIM. When a symmetrical pulse passes through 
a NZPIM slab near the zero-refractive-index point, its spatio- 
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temporal spectrum will be filtered asymmetrically for negative 
and positive refractive index regions, which results in the ZB- 
like effect in time domain. The strength of tailed oscillations 
depends on the sample thickness h, temporal half-width T, and 
the spatial half-width Wo while the oscillating frequency only 
relates to T and Wo, and do not change with h. The results 
we have got are very similar to what the authors of Refs. ]8[] 
and Hill obtained based on the concepts of relativistic quan- 
tum mechanics, however, our analysis is a classical method 
which is irrelevant to the relativistic quantum mechanics. 
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